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The low temperature magnetization process of antiferromagnetic spin-S 1 chains doped with mobile 
spin-(5 — 1/2) carriers is studied in an exactly solvable model. For sufficiently high magnetic fields 
the system is in a metallic phase with a finite gap for magnetic excitations. In this phase which 
exists for a large range of carrier concentrations x the zero temperature magnetization is determined 
by x alone. This leads to plateaus in the magnetization curve at a tunable fraction of the saturation 
magnetization. The critical behaviour at the edges of these plateaus is studied in detail. 
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Synthetization of new magnetic materials and avail- 
ability of very high magnetic fields provide new pos- 
sibilities to study the magnetization process of low- 
dimensional quantum spin systems. In particular so- 
called spin liquids realized in quasi-one dimensional an- 
tiferromagnetic systems such as spin chains, spin ladders 
and exchange-alternating spin chains attract much inter- 
est at present due to the possible occurence of magne- 
tization plateaus associated with gapped excitations. In 
addition to saturated magnetization M s such plateaus, 
i.e. regions where the magnetization does not depend 
on the magnetic field for sufficiently low temperatures, 
are admissible from topological considerations at certain 
fractions of M s depending on the value of the spin of the 
substance and the translational symmetry of the ground 
state. Necessary conditions for the occurence of plateaus 
have been formulated by Oshikawa et al. jl| employing 
a generalization of the Lieb-Schultz-Mattis theorem: for 
a spin-5 chain with a magnetic unit cell containing q 
magnetic moments this feature can appear at rational 
values (M) with integer q(S — (M)). The existence of 
these phenomena in a variety of models has been estab- 
lished by numerical and analytical studies of various low- 
dimensional magnetic insulators including spin chains, 
spin ladders and systems with multi spin exchange or 
exchange anisotropies [pl-jiUf. Very recently, several ex- 
perimental observations of such magnetization plateaus 
at non-zero (M) have been reported 

A common feature in these systems is that the plateaus 
in the magnetization curves appear at certain simple frac- 
tions of the maximal value M s as a consequence of their 
topological origin. In this letter we report on a mecha- 
nism leading to gaps for magnetic excitations at magneti- 
zations which can be controlled by suitable preparation of 
the sample, namely doping. We study this phenomenon 
in the framework of a recently introduced class of inte- 
grable models for doped Heisenberg chains which may 
be used as a basis for studies of certain features of doped 
transition metal oxides 17 lq]. Starting from the double- 



exchange model 19 1 a strong ferromagnetic Hund's rule 



and localized quantum spins (5—1/2) arising from the ti g 
electrons allows to introduce an effective Hamiltonian on 
a restricted Hilbert space with maximally allowed spin S' 
on a given lattice site (2^j2l[], i.e. S' = S if the electronic 
state on this site is occupied, or S' — S — 1/2 if there is 
no e g electron (denoted as a hole in the following). This 
derivation of a low energy Hamiltonian generalizes that 
of the t-J model from the Hubbard model (22) which 
is corresponds to the case of S = 1/2, i.e. no localized 
spins. Numerical studies of the S > 1 variants of these 
models have been performed to gain a better understand- 
ing of experimental findings for the doped Haldane sys- 
tem Y2- 2; Ca 2 ;BaNi05 (S — 1) [^0) and manganese oxides 
such as Lai^Ca^MnOa (S = 2) ||. 

Below we consider integrable models of this type in one 
spatial dimension. Similar to the models obtained from 
the general procedure outlined above their Hamiltonians 
are of the form 
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Here X^p and T^' describe the (antiferromagnetic) ex- 
change and hopping of the holes between sites i and j of 
the lattice, respectively. SU(2) invariance of the model 
implies that they can be written as polynomials of the op- 
erator Si-Sj (see e.g. ^T|). The form of these polynomials 
is fixed in the integrable models JT8fl . For example, in the 
5=1 case with possible relevance to the doped Nickel 
oxides the antiferromagnetic exchange terms are given in 
terms of bilinear and biquadratic Heisenberg couplings 
depending on the values Sij 6 {1/2, 1} of the spins on 
sites i and j pjv : 
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(Note that the undoped chain, i.e. Si = 1 for all i, is 



coupling between the spins of the itinerant e g electrons 



the integrable spin-1 Takhtajan-Babujian model [£4 25 
while the completely doped chain is the spin- 1/2 Heisen- 
berg chain with bilinear exchange.) Similarly, the kinetic 
term of the integrable spin-1 model reads 
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where Vij is an operator permuting the states on sites i 
and j thereby allowing the spin-1/2 "holes" to propagate. 
The additional exchange term in this expression leads to 
different hopping amplitudes t(Sij) depending on the to- 
tal spin Sij on the participating sites, i.e. £(1/2) = — 1 
and £(3/2) = +1/2 for the possible values Sij — 1/2 and 
3/2, respectively. These amplitudes differ from the values 
proposed in Ref. [^0| for the doped Nickel oxide, namely 
£(1/2) = -1/2, £(3/2) = 1. This is one reason for the 
absence of a ferromagnetic phase in the integrable model 
(U) (sec JltJ for a discussion of the other differences). 

For general S the integrable models are constructed 
from solutions of a Yang-Baxter equation and can be 
solved by means of the algebraic Bethe Ansatz [|l8| . Their 
thermodynamical properties at finite temperature T can 
be obtained from the solution of the thermodynamic 
Bethe Ansatz (TBA) equations, i.e. the following set of 
coupled nonlinear integral equations 

e„(0 = Ts * ln[l + e e - 1 «)/ T ][l + e ^+^ T ] 
-2nS 1h2S s(0 - 6 nA Ts * ln[l + e -^/ T ] , 

(2) 

-[2ira 2S * s(£) + /i] - Ts * ln[l + e eiK)/T ] 
= k(0 + TR* ln[l + e~ K( « )/T ] . 

Here (/ * g) (£) denotes a convolution in the space 
of rapidities £, a„(£) = (2n/ir) (4£ 2 + n 2 )~\ s(£) = 
(2 cosh7r£) _1 and R — a 2 * (1 + (I2) -1 . Eqs. (^) are to 
be solved subject to the condition lim„_ +00 (e„/n) = H 
with the external magnetic field H and \i is the chemical 
potential for the holes controlling their concentration. In 
terms of the functions e n (£) and k(£) the free energy of 
this system reads (£q is the ground state energy of the 
spin-S Takhtajan-Babujian chain for H = |pq| ) 

~F(T, H,ri = ^--T j s(0 ln[l + e ea *<«/ r ] 

-T | d£ (a 2S * a) (0 Ml + e" K («)/ T ] . (3) 

The low temperature (iJ 3> T) phase diagram for the 
spin-1 system has been obtained in Ref. []l7| , qualita- 
tively the same behaviour is found for general S > 1 [18| . 

Here we study the properties of these systems in a 
magnetic field at fixed doping. For hole concentrations 
< x < x c (S) (see Fig. [l]) the low energy excitation 
spectrum of the system allows to identify four interme- 
diate field phases (labelled A, Bi, C and B2 in Fig. 1 
of Ref. |Q ) for < H < H s before the system is com- 
pletely polarized for H > H s . Of particular interest in 
the present context is the phase C: since the system is not 
ferromagnetically polarized one expects nontrivial exci- 
tations for both the charge and magnetic degrees of free- 
dom. In the neighbouring phases Bi 2 these excitations 



are massless leading to an effective description of these 
phases in terms of a Tomonaga-Luttinger model. The 
analysis of the T = limit of the TBA eqs. (||) shows 
that in phase C only one of these modes is gapless Q . 
The resulting low-energy theory is that of a single mode 
with dispersion 

e2s(0 + f Q d£> K(£ - £)e 2S (?) = 4s (0 (4) 

J-Q 

where K(£) = 2£ 2 ±7 1 a 2fc (£), 4s(0 = {2S - H - 
[X — 2Tr^2 k=1 a 2 k-i(l;) and Q is a function of mag- 
netic field and chemical potential through the condi- 
tion £2s(±Q) = 0. The corresponding hole concentra- 
tion x — J_q d£ <72s(£) is obtained from an equation 
for <7 2 s similar to §Q) with driving term replaced by 

°2s(0 = Sfc=i a 2fc-i(0- Further analysis of the zero 
temperature limit of the TBA eqs. (||) shows that the 
massless mode in this phase carries the charge degrees 
of freedon while all magnetic excitations are gapped, i.e. 
k(£) < and e„^2s(£) > for all £. Hence, x and the 
magnetization M p = S — 3x/2 are constant throughout 
this phase for fixed Q, i.e. (2S — 5) H — fi = const. This 
implies plateaus in the magnetization curve M(H) below 
the saturated value M s = S — x/2 (see Fig. || for S = 1). 
The end points of these plateaus are H c \ = —2\i and 
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As H — > H c \ 2 from inside the plateau region the spin 
gap closes as A oc \H — H c i 2 \. 

For finite temperatures the full set of TBA eqs. has 
to be solved to determine the magnetization curves. In a 
sufficiently strong magnetic field H 3> T however, the en- 
ergies e„>2S are gapped and can be eliminated from the 
TBA eqs. (j|) @. For the doped S = 1 chain this pro- 
cedure leads to a coupled set of three nonlinear integral 
equations which are straightforward to solve by iteration. 
Choosing the chemical potential such that the hole con- 
centration x = —dF/dfi is fixed the magnetization and 
magnetic susceptibility can be obtained from the ther- 
modynamical potential f2(T, H, x) = F(T, H, /1) + \xx. In 
Fig. |3| we present the resulting data for various tempera- 
tures. They clearly show the formation of plateaus with 
decreasing temperature and the singular behaviour aris- 
ing in the vicinity of the transitions into the spin gap 
phase as expected from the analysis of the zero temper- 
ature phase diagram. 

Remarkably, the nature of these singularities on the 
two critical end points H c i 2 of the plateau is quite dif- 
ferent: for the magnetic insulators discussed in the in- 
troduction the singular part of the magnetization near 
the plateaus has been predicted to show a square root 
behaviour M due to the similarity of the transition to 
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a commensurate-incommensurate transition p6j . A reli- 
able numerical verification of this prediction — even for 
an integrable model — is extremely difficult for transi- 
tions other than the one into the ferromagnetically po- 
larized state JjJ. In the model considered here such dif- 
fulties arise for H — > H c i only: for H < H c i the zero 
temperature magnetization shows a critical behaviour 
oc (H c i — H) a consistent with the square root behaviour 
a = 1/2 within the numerical accuracy of our data. On 
the other hand, near H C 2 the magnetization depends lin- 
early on the external field, i.e. M — M p cx H — H C 2 for 
H > H C 2. This difference in the critical behaviour is also 
evident in the temperature dependence of the magnetic 
susceptibility near -ff c i.2, in particular \ ss const, at the 
high- field end of the plateau (see Fig. ||(b)). Note, that a 
similar T-dependence has been observed in experiments 
on certain spin- 1/2 Heisenberg ladders |2^j2^] . 

This different singular behaviour is a consequence of 
the coupling between the two massless excitations present 
in the Tomonaga-Luttinger phases outside the interval 
H c i < H < H C 2'. without the magnetic field breaking 
the spin-5f7(2) these modes can be assigned usually to 
spin and charge excitations separately based on their dif- 
ferent symmetries. In an external field, however, this 
assignment leads to a coupling of the two sectors. In cer- 
tain cases this interaction can be removed by allowing 
for mixing of the corresponding quantum numbers [^9| . 
Analysing the Bethe Ansatz wave functions we can de- 
termine the relation between the total charges Qi,2 in 
the two gapless sectors to the physical quantum num- 
bers, i.e. number of holes and z-component of the total 
spin. In the phase for H < H c \ a change in magneti- 
zation at fixed doping affects the total charge in one of 
the sectors only. Necglecting the coupling to the other 
sector, bosonization then gives the familiar square root 
singularity of the magnetization at H < H c \ p6|| . For 
H > H C 2, however, a change of magnetization requires 
AQi — — AQ2 for fixed doping x. As the critical field 
H C 2 is approached from above, this fact together with 
the coupling between the two sectors leads to a disper- 
sion ei(fc) w (vpk 2 /4A) — A of the "incommensurate" 
soft mode with A cx (H — H C 2) 2 rather than the usual 
behaviour A oc \H — H C 2\ for fermions. This field depen- 
dence immediately gives the linear 77-dependence of the 
magnetization for H > H C 2- 

In summary we have studied the properties of doped 
Heisenberg chains in a magnetic field in the framework of 
a Bethe Ansatz solvable model. We find plateaus in the 
magnetization curves at certain values M p of the mag- 
netization. To our knowledge this is a novel feature in 
an integrable model thereby providing the basis for more 
detailed studies aiming at a better understanding of the 
mechanism for the occurence of these plateaus and the 
critical behaviour at their end points. M p can be contin- 
uously tuned by changing the concentration of carriers. 



This is a fundamentally different from the fixed values 
of M p obtained from topological arguments in the mag- 
netic insulators studied previously. An extension of such 
arguments to the case of doped chains might be possi- 
ble within a classical treatment of the localized t2 g spins 
in the double-exchange model: in this limit the ground 
state has an incommensurate magnetic structure with kp 
periodicity ]|(|. However, unlike the models considered 
here this approach also leads to gapped charge excita- 
tions. An alternative attempt to describe the plateaus at 
tunable fractions of the maximal magnetization within a 
bosonized model relies just on this existence of a gapless 
charge mode (3l]]. Similarly, the second feature where 
the plateaus discussed here differ from the ones in quan- 
tum spin chains and ladders can be understood only as a 
consequence of the existence of a second massless mode: 
the relation of the physical quantum numbers to the con- 
served charges of the effective low energy theory together 
with the coupling of the gapless channels conspire to give 
the linear field dependence of the singular part of the 
magnetization observed near the critical point H c i- Note 
that this mechanism does not restrict the critical expo- 
nents to the values 1/2 and 1 discussed in this letter. 

Further studies of the low energy properties — in par- 
ticular the analysis of the asymptotics of correlation func- 
tions as H c \ ^2 are approached — in this solvable micro- 
scopic model will lead to new insights into the critical 
behaviour at the plateau transitions and possible the re- 
lated ones into Mott insulating phases of interacting par- 
ticles. Furthermore, the phenomena reported in this let- 
ter may be verified in experimental studies of the magne- 
tization process in the doped, effectively one-dimensional 
transition metal oxides mentioned above. 
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FIG. 1. Maximal concentration of holes x c (S) for the exis- 
tence of a magnetization plateau vs. 5* -1 (the line is a guide 
to the eye only). 
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FIG. 2. Zero temperature magnetization curve of the 
doped 5=1 chain for the hole concentrations x = 0.1, 0.2, 
0.3 (top to bottom). 
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FIG. 3. Magnetization curve (a) and magnetic susceptibil- 
ity (b) of the doped S = 1 chain with hole concentration 
x = 0.2 at temperatures T = 0.02, 0.05, 0.1, 0.2 (in units of 
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